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approximation the tensor amplitude for longitudinal polarization of the virtual photon is finite, 
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divergence has zero projection on the polarization vector of the final photon, so that the physical 
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I. INTRODUCTION 



Deeply virtual Compton scattering (DVCS) is currently receiving a lot of attention, as an exclusive process 
closely related to inclusive deep-inelastic scattering (DIS) The theory of DIS is well understood and can be 

formulated both in the field-theoretical framework of QCD or the particle-based language of the parton model. 
The relevant information about the non-perturbative structure of the nucleon is contained in expectation values 
of QCD light-ray operators in the nucleon, which can be interpreted as the parton distributions. Similarly, the 
DVCS amplitude in the asymptotic region can be expressed through off-diagonal (non-forward) matrix elements of 
light-ray operators, which also possess a partonic interpretation in terms of so-called "skewed" parton distributions 
(SPD's). The initial motivation for considering DVCS was the hope that by measuring the SPD's one could obtain 
information about the angular momentum distributions of quarks and gluons in the nucleon [jjj. However, it 
was soon realized that the DVCS process is of interest in itself, independently of its connection to the nucleon 
spin problem. It is the simplest representative of a whole class of exclusive electroproduction processes, including 
also hard meson production, for which factorization can be proven in QCD at leading twist level [0j3]. The 
first experimental observations of DVCS have recently been reported by the ZEUS ||, HI and HERMES || 
collaborations, and an extensive experimental program is planned at JLAB H. The DVCS process is also related, 
by crossing relations to the electroproduction of pairs of hadrons off a real photon, such as 7*7 — > 2-zr fl0| , 

which can be studied at e + e~ colliders. 

The virtual Compton amplitude in the asymptotic region of large virtuality of the incoming photon, — qf — > 00, 
and small momentum transfer squared, t <C 1 GeV 2 , has been studied in QCD using expansion techniques familiar 
from DIS: The collinear expansion in momentum space p"l) , and the non-local light -cone expansion in coordinate 
space Jl^ , 0|. Originally only the contribution from twist-2 operators was included in the calculation of the tensor 
amplitude [Q,^) . This approximation does not include subleading terms proportional to the transverse part of the 
momentum transfer, rj_ = A, and the tensor amplitude thus obtained formally violates electromagnetic gauge 
invariance (trans versality) by terms linear in To remedy this, Guichon and Vanderhaeghen Jl5| proposed to 
add a term linear in A, producing a tensor amplitude which is transverse up to terms of order A 2 ~ t. It is clear 
that in a QCD treatment such a term could come only from contributions of operators of higher twist (> 2) . Indeed, 
calculations by various groups |l6|-|ll| have recently shown that the term introduced in Ref. jl5| appears naturally 
as part of the twist-3 contribution to the Compton amplitude. Anikin et al. Pq] , applying the momentum-space 
collinear expansion up to twist-3 level, reproduced not only the "improved" twist-2 structure of Ref. fLq) , but 
found also other, separately transverse twist-3 structures in the tensor amplitude, which can contribute to the 
observables of the DVCS process. A similar result was obtained by Penttinen et al. |T7| , who computed the DVCS 
amplitude to accuracy 1/ \/ —q\ in a parton model approach. Within the framework of the coordinate-space light- 
cone expansion the twist-3 contribution to the DVCS amplitude was studied by Belitsky and Miiller In 
particular, these authors demonstrated that in order to get a gauge invariant result up to terms of order tj q\ it 
is sufficient to retain only a certain part of the full twist-3 SPD's which is obtained by Wandzura-Wilczek-type 
relations from the twist-2 distributions. Kivel et al. |2(J noted that the Wandzura-Wilczek expression for the 
tensor amplitude contains a divergence in the part corresponding to transverse polarization of the virtual photon. 

In this paper we address the question of transversality of the DVCS amplitude within the framework of the 
light-cone expansion of the current product in QCD string operators in coordinate space developed by Balitsky 
and Braun jl3|]. (A short summary of our results has already been given in Ref. |L9]].) The investigation consists 
of two main parts. First, we show that transversality of the light-cone expansion can be achieved, in the minimal 
sense, by including in the expansion certain operators which are algebraically of twist 3 but given by total derivatives 
of twist-2 string operators. We shall refer to these as "kinematical" twist-3 contributions, in order to distinguish 



*BKimlein and Robaschik jl4| studied instead of the tensor DVCS amplitude the individual photon-hadron helicity ampli- 
tudes. The problem with the gauge invariance of the twist-2 contribution does not arise as long as one considers only the 
leading helicity amplitudes in 1/ yj— qf. 
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them from the twist-3 contributions involving quark-gluon operators which cannot be reduced to total derivatives 
by the equations of motion. The operations necessary to include the "kinematical" twist-3 contributions in the 
light-cone expansion can be performed using the calculus of QCD string operators. The result has the same form 
as the usual twist-2 part of expansion, but with certain finite translations of the center coordinates of the twist-2 
string operators (similar expressions were derived independently in Refs. Jl^j20[| ). 

Second, we investigate whether the approximation of keeping only the "kinematical" twist-3 contributions re- 
quired by transversality leads to useful results at the level of the hadronic DVCS amplitude. To this end, we 
parametrize the matrix elements of the fundamental twist-2 operators in terms of two-variable spectral functions 
(double distributions) ||f21 22 j|. The Compton amplitude in the DVCS limit, including the "kinematical" twist-3 
contributions, can be expressed in terms a number of one- variable skewed distributions (SPD's), which are obtained 
in the form of one-dimensional reductions of the twist-2 double distribution. The resulting tensor amplitude is 
manifestly transverse up to terms of order tjq\ or m 2 /q 2 (m is the target mass), i.e., its contraction with the 
four-momentum of the initial photon, qi^T^, and that of the final photon, T^ u q2 U , are both of order tjq\ or 
m 2 /qf. We find that our approximation gives a finite result for the part of the tensor amplitude corresponding 
to longitudinal polarization of the initial photon, while the part corresponding to transverse polarization contains 
a divergence. However, as we shall show, the divergent part has zero projection on the polarization vector of the 
final (real) photon, so it does not affect the physical helicity amplitudes. 

Our approach of keeping only the "kinematical" twist-3 operators in the coordinate-space light cone expansion 
is actually equivalent to the Wandzura-Wilczek-type approximation for the twist-3 SPD's introduced in Ref. |l8| . 
In addition to giving an alternative derivation of this approximation, an important advantage of our treatment 
is the possibility to incorporate the formalism of double distributions ]^, ^l] , p2] ,p[ . The twist-2 as well as the 
new twist-3 SPD's are obtained by reduction formulas from the fundamental twist-2 double distribution. This 
allows for a straightforward derivation of the properties of these new functions and the Wandzura-Wilczek-type 
relations between them. In particular, the fact that this approximation gives finite results for the tensor amplitude 
for longitudinally polarized initial photon, but contains a divergence in the case of transverse polarization (which, 
however, drops out of the physical helicity amplitudes), can immediately be understood on the basis of the reduction 
formulas. 

The plan of this paper is as follows. Section || gives a brief summary of the kinematics of virtual Compton 
scattering in the generalized Bjorken limit. In Section III we discuss the peculiarities of the light-cone expansion 
in the "non-forward" case, where total derivative operators contribute. Our treatment is based on the formalism 
developed by Balitsky and Braun [jl3| . In Section [II A we recall the main steps in the twist decomposition of string 
operators. In Section IIIB we investigate the problem of transversality of the light -cone expansion at operator 
level in coordinate space. We observe that the twist-2 contribution is not transverse, and show that transversality 
can be restored by including total derivatives of twist-2 operators, which are algebraically of twist 3. We first 
formulate this as as an iterative procedure generating an infinite series of total derivatives of twist-2 operators, 
then show that this series can be summed up in closed form to give finite translation operators acting on the twist-2 
string operators (details are given in Appendix |A|) . A manifestly transverse expression for the light-cone expansion 
of the current product in terms of twist-2 operators and their total translations is presented. In Section [V we 



use our result for the light-cone expansion including "kinematical" twist-3 terms to compute the hadronic DVCS 
amplitude. For simplicity we consider the case of a pion target, which has spin and definite C-parity. In 
Section IV A we describe the spectral representation of the fundamental twist-2 matrix elements in terms of double 
distributions. The matrix elements of the vector-type string operators including the "kinematical" twist-3 terms 
are then derived in Section IV B. The new SPD's are introduced through reduction formulas, and their properties 
are discussed. The DVCS amplitude for the pion is computed in Section IV C. We inspect the singular integrals 
resulting from the inclusion of the "kinematical" twist-3 terms. We show that the tensor amplitude for transverse 
polarization of the initial (virtual) photon contains a divergence, but that this divergence is "harmless" in the 
sense that it has zero projection on the polarization vector of the final photon. In Section IV D we reformulate 
our results for the matrix elements and the Compton amplitude, introducing a new SPD as a Wandzura-Wilczek- 
type transformation of the basic twist-2 SPD. This allows us to demonstrate that our approach is equivalent 
to the Wandzura-Wilczek approximation proposed in Ref. |lq j. In this language, the singularity for transverse 
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polarization described in Section IV C| appears due to a discontinuity of the Wandzura-Wilczek-transformed SPD 
at x = £, which occurs independently of the dynamical behavior of the underlying twist-2 matrix elements. Our 
conclusions are summarized in Section [v|. 

In this paper we study the properties of the DVCS amplitude including "kinematical" twist-3 contributions 
at tree level. We shall not be concerned with logarithmic corrections resulting from the scale dependence of the 
operators (the evolution of the skewed distributions). For the twist-2 operators this problem has exhaustively 
been treated in the literature, see e.g. Refs. p|,p|,p3-27|. The generalization of the evolution equations to the 
"kinematical" twist-3 contribution discussed in the present paper is in principle straightforward but will not be 
pursued here. 



II. VIRTUAL COMPTON AMPLITUDE IN THE GENERALIZED BJORKEN LIMIT 

Virtual Compton scattering. For a unified description of the kinematics of processes such as DIS and DVCS let 
us imagine a general virtual Compton scattering process off a hadron, 

7*(<7i) + HPi) - 7*fe) + h(p 2 ), (2.1) 

where the incoming photon has space-like virtuality, q 2 < 0, and the final photon may have any virtuality allowed by 
four-momentum conservation; in particular, we shall allow for the final photon to be real, q 2 = 0. It is convenient 
to use as independent momentum variables the average of the photon and hadron momenta, q and p, and the 
momentum transfer, r: 

r r 
Pi,2 = P±o' 91,2 = qT 2' *- 2 ' 2 ^ 

From the mass shell conditions, p\ = p\ = to 2 , it follows that 

(pr) = 0, p 2 = m 2 - -, t = r 2 . (2.3) 

In the following we shall neglect the target mass and put to = 0. 

The amplitude for the virtual Compton process is defined by the transition matrix element of the time-ordered 
product of two electromagnetic current operators between the hadronic states: 

= i J d A ze^(p-r/2\TJ^-z/2)J v {z/2)\p + r/2). (2.4) 

This tensor is transverse with respect to the incoming and outgoing photon momenta: 



r 

q -2 



v = o, T ^{ q+r i)u = °' (2,5) 



which is an immediate consequence of current conservation, d^J^x) = 0{j] Separating the symmetric and antisym- 
metric tensor parts of T M „ (with otherwise same arguments), 



-(T ltu ±T vli ), (2.6) 



the two conditions Eq.(2.5) can also be stated in the form 



' Strictly speaking this is correct only up to terms arising from the differentiation of the step functions in time occurring in 
the reduction of the time-ordered product ("seagull terms"). We shall not consider these contributions, since they do not 
play a role at large q 2 . 
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% T [nv] = y T {^}- ( 2 - 8 ) 

One notes that in the case of non-forward scattering, r ^ 0, the transversality conditions generally imply relations 
between the symmetric and antisymmetric parts of the amplitude. It is only in the limit of forward scattering, 
r = 0, (as is relevant e.g. to deep-inelastic scattering) that the two transversality conditions decouple and turn 
into conditions on the symmetric and the antisymmetric parts separately. 

Generalized Bjorken region: DIS and DVCS. We shall consider the amplitude for the virtual Compton process, 
Eq.( [2.l| ), in two different asymptotic regions, both characterized by large virtuality of the incoming photon. One 
is the region corresponding to deep-inelastic scattering (DIS), where one is interested in the forward scattering 
amplitude (q± = q2 = q,Pi = Vi — P> r — 0) in the limit q\ — > oo, (pi<7i) — * oo, with xb = —qf/[2(piqi)] fixed. 
The other is the region of deeply virtual Compton scattering (DVCS), where one is dealing with the amplitude 
for production of a real photon, q\ = 0, in the limit qf — > oo, (piqi) — > oo, again with xb = — <Zi/[2(pi?i)] fixed. 
This process requires a non-zero momentum transfer from the initial to the final hadron, (pi — p2 = r ^ 0), with 
(rqi) — > oo proportional to q 2 . One considers a situation where t = r 2 < does not grow with q\. In order to treat 
both asymptotic regions within a common framework it is convenient to introduce two scaling variables, 

where rj is referred to as "skewedness" . The virtualities of the incoming and outgoing photon can then be expressed 
as 



qf= [l + jjq' + ~, (2.10) 



|1-|)? 2 + ^ (2.11) 



and the Bjorken variable becomes 



x b = W7- : \ = ttz.- ( 2 - 12 ) 



2{pxqi) I + 77' 



Because of Eq.(2.1C) we may generally use q 2 instead of q\ as the large invariant. The two situations, DIS and 



DVCS, now correspond to^] 

DIS: x B = f, 



2£ (2-13) 

DVCS: V = £, xb — 



Aim of the present calculation. In this paper we calculate the virtual Compton amplitude, Eq.(2.4), in DVCS 
kinematics, neglecting contributions of order t/q 2 or m 2 j q 2 . Specifically, we aim to compute the tensor amplitude 



to an accuracy such that the transversality conditions, Eqs.(2.5) viz. (2.8), are satisfied up to terms of order t/q 2 
or m 2 j q 2 . In this approximation the kinematical limits for £ are 

0<£<1, (2.14) 

which is equivalent to < xb < 1. For DVCS it is convenient to separate from the momentum transfer, r, the 
component parallel to p, which is responsible for producing the large scalar when contracting with q, i.e., to write 



Because of the identification of £ with n in DVCS kinematics, £ in DVCS is frequently also referred to as "skewedness" . 



G 



r = 2£,p + A. 



(2.15) 



The remaining part of the momentum transfer, A, satisfies (Aq) = — 1/4 and (Ap) = — £t/2. Since p is time-like 
and q space-like, and since furthermore A itself is space-like, A 2 = (1 — £ 2 )i < 0, it follows that all components of 
A must be of order \t\ x / 2 : 

A cc \t\ x ' 2 . (2.16) 

This means that when computing the DVCS amplitude using the light-cone expansion we may expand the ingredi- 
ents — the momentum-space quark propagators and the hadronic matrix elements — in A, since this four-vector 
cannot produce any large scalars even when contracted with the photon momentum, q. This would not be possible 
for the full momentum transfer, r, since (rq) is not small. In the language of the parton model, the component 2£p 
of the momentum transfer would be identified with the longitudinal component, while A would be the transverse 



component, r±. Note that in our approach the decomposition of the momentum transfer, Eq.(2.15), appears on 
purely kinematical grounds, without reference to a particular frame. 



III. LIGHT-CONE EXPANSION IN THE NON-FORWARD CASE 
A. Twist decomposition of string operators: Twist 2 



The behavior of the virtual Compton amplitude, Eq.( [2.4|) , in the generalized Bjorken limit (DIS or DVCS) is 
dominated by the light-cone singularity of the time-ordered product of the two electromagnetic currents, which is 
a consequence of the asymptotic freedom of QCD. An asymptotic expansion of the Compton amplitude in inverse 
powers of 1/q 2 can be derived from the light-cone expansion of the current product. In this section we discuss the 
peculiarities of the light-cone expansion in the "non- forward" case, i.e., finite momentum transfer in the matrix 
elements, as is relevant e.g. for DVCS. An additional complication compared to the "forward" case arises due to 
the fact that total derivatives of operators can contribute to the matrix elements, which greatly enlarges the set 
of operators participating in the expansion to a given accuracy. In particular, we show here that keeping the total 
derivatives of twist-2 operators (which are operators formally of twist 3) is essential for maintaining transversality 
of the leading term in the light-cone expansion. We shall proceed step by step. We first investigate the well-known 
twist-2 contribution and show that, in the non-forward case, it is not transverse to the accuracy required. We then 
formulate a procedure how transversality can be "repaired" by including total derivatives of twist-2 operators, 
order-by-order in the number of total derivatives. Finally, we show that the series of total derivatives of twist-2 
operators can be summed up in closed form, resulting in an expression involving certain finite translations of the 
original twist-2 operators, which is manifestly transverse. 

Expansion in string operators. We shall be concerned with the light-cone expansion of the current product, 

Knv(x,y) = iTJ ti (x)J v (y), (3.1) 

where J^ix) is the vector current operator. It is convenient to introduce the center and relative coordinates of the 
two points, x and y, 

X = (x + y)/2, z = y-x, = Itz/2, (3.2) 

and to regard the current product as a function of these variables. The light-cone expansion is derived from a 
formal expansion of the current product in QCD string (non-local) operators in a background gauge field, which are 
subsequently decomposed in operators of definite twist. This expansion is generated by contracting the quark fields 



in Eq.(B.l) with the quark Green function in a background gauge field, and expanding in powers of insertions of the 
background field in the Fock-Schwinger gauge. The resulting expansion can be formulated in a gauge-invariant way 
by introducing QCD string operators with fields connected by gauge links. It is known that the leading light-cone 
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singularity is contained in the contraction corresponding to the simplest two-point string operator (the "handbag" 
diagram) : 

U^(z\X) = iTJ^X - z/2)J v (X + z/2) 

= ${X - z/2)g[l - z/2, X + z/2]ij(X + z/2) 
+ [z — > —z, /J, v). 

+ terms 1/z 2 ,... (3.3) 



Here 



2ir 2 (z 2 -i0) 2 J (2tt) 4 k 2 + iO 



2 = ^ fifc . ) (3 . 4) 



is the free quark propagator in coordinate representation. Reducing the product of three gamma matrices, Eq.(3.3) 
can be written as 



n^{z\x) 

2z 



where 



ma [i>{X - z/2) la [X - z/2,X + z/2]iP{X + z/2) - {z -> -z)] 

pva [iP(X - z/2)w 5 [X - z/2,X + z/2]iP(X + z/2) + ( z -> -z)]} , (3.5) 



Supvo = 7*r [j^jpjuja] = g^pQva - 9 ,i, v g pa + 9p,a9vp- (3.6) 



In the string operators in Eq.(3.3) and ( |3.5| ) the gauge link runs along the straight line connecting the two space- 
time points: 

[x,y] = Pexpz / du (x ~ y) a A a (ux + uy), u = l — u. (3-7) 



For brevity we shall in the following not write the links between the fields; it will always be understood that the 
fields in the string operator (including possible insertions of gauge fields at intermediate points) are connected by 
gauge links along the straight line. 

Differentiating string operators. For a given background gauge field the link operator defined by Eq.(3.7) is a 
uniquely defined function of the end points. In particular, when differentiating the string operator with respect 
to its end points one is changing the contour of integration by an infinitesimal amount, so the link operator gets 
differentiated as well. There are generally two contributions from differentiating the link operators. One gives the 
gauge potentials at the end points, which combine with the derivatives acting on the quark fields to give covariant 
derivatives. The other can be expressed as a line integral of the gauge field along the contour of the string. For 
the derivative of the string operators in Eq.(3.5) with respect to the relative coordinate we find: 



~ $(X - z/2) (v„ la + la V„) il>(X + z/2) 



1 

dtt 4>{X - z/2) la z a F^ a (X + tz/2)ip(X + z/2), (3.8) 

-l 

where gauge links between the fields are implicit everywhere; for the string operator with Dirac matrix 7^75 one 
should replace la ~> 7o-75 everywhere. The derivative with respect to the center coordinate (we refer to it as 
"total" derivative) takes the form 



8 



ax, 



4{X - z/2)i^{X + z/2) 



<P(X - z/2) (- V M la + la V M J r/>(X + z/2) 
2/ / dt ^{X - z/2)iaZ a F tia (X + tz/2)iP{X + z/2). 



Here 



V M i>(x) = dnip(x) - iA fl (x)ip(x), 
4>(x) V M = —dn%j}{x) - i'4>{x)A ll {x) 

are the right and left covariant derivatives^, and 



F u 



i Vu,Vq 



(3.9) 

(3.10) 
(3.11) 

(3.12) 



the field strength. 

Twist-2 part of string operators. The string operators in Eq.( |3.5[ ) have as yet no definite twist. Consequently, 
the "handbag" contribution to the current product, Eq.(^), as it stands, contains also subleading contributions in 
1/z 2 , which would be of the same order (l/z 2 or less singular) as contributions from quark-gluon string operators 
not included in Eq.( |3.5| ). The leading l/z 4 singularity is given by the twist-2 part of the string operators. The 
latter is defined as those operators obtained by formally Taylor-expanding the string operators in Eq. ( |3.5| ) in the 
separation, z, and retaining only the totally symmetric traceless parts of the coefficients in the expansion:Fj 



$(X-z/2)ia^(X + z/2)] twist 2 

1 

= — z ai ■ ■ ■ z«„ ${X) i{a Voi • ■ ■ V a „) - traces ip(X), 



n=0 



(3.13) 



and similarly for the operator with Dirac matrix 70-75. Here V=(V + V)/2. As was shown in Ref. [Q, the 
two operations, "symmetrization" and "subtraction of traces" , can in fact be carried out directly at the level of 
non-local operators, without having to go through the series of local operators p8j . The part of the string operator 
corresponding to totally symmetric local tensor operators is projected out by the operation 



$(X-z/2)iaiP(X + z/2)Y 



d f 1 - 
— / dt ip(X -tz/2)ziJ>(X + tz/2). 

<JZ a Jg 



(3.14) 



The subtraction of traces can be achieved by noting that the tracelessness of the local operators in the expansion, 
Eq.( [3.13| ), implies that the string operator contracted with z a should satisfy the d'Alembert equation with respect 
to z: 



(X - z/2)zip(X + z/2)] 



traceless 



0. 



(3.15) 



Thus, combining the two operations, the twist-2 part of the vector string operator (for arbitrary z 2 ^ 0) is obtained 
by substituting in the R.H.S. of Eq.(3.14) the solution of Eq.( |3.15 ) with "initial conditions" given by the operator 
ijj(X — z/2)zijj(X + z/2) at z 2 = 0. On the light cone (z 2 = 0) trace subtraction is irrelevant, and the twist-2 part 



of the vector string operator coincides with the symmetric part as defined by Eq.(3.14). 



" Since we are dealing with derivatives of fields of modified arguments it is important to maintain the distinction between 
the gradient evaluated at the point z/2, d fl tf)(z/2) = (d/dxi_ l )ifj(x) x=z /2, and the derivative (d/de fl )ip(z/2) = (1 / '2) d^ij) {z/2). 

"Note that this definition is at a purely algebraic level and does not imply any statement about the convergence of the 
series of local operators. 
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When parametrizing the hadronic matrix elements of the operators, the z 2 - ("trace") terms in the matrix 
element are proportional to dimensionful scalars characterizing the target — the target mass, m 2 , or the square of 
the momentum transfer, t. The inclusion of the z 2 terms in the calculation of the Compton amplitude then generates 
m 2 /q 2 - (target mass) and £/g 2 -corrections to the DVCS amplitude. Here we shall neglect such corrections, so that 
the z 2 -terms in the matrix elements need not be subtracted explicitly. 

Twist-2 contribution to the light-cone expansion. The twist-2 part of the current product, which we denote by 
11^2, is given by Eq.( |3.5[ ), with the string operators replaced by their symmetric parts (we do not worry about 
trace subtraction): 



2z n 



3 [LpVO 



(X-z/2) la ^(X + z/2) -{ z ^-z)Y 



-kppu, [i>(X ~ z/2) laTo ^(X + z/2) + (z -f -z)] y } 



(3.16) 



Using the property Eq.( 3.15 ), and the fact that the coefficient function — the free fermion propagator, Eq.(3.4) 
satisfies 



one can easily show that the twist-2 part of the current product satisfies 

—n m - o 

where 



—Tl (0) 



n 



(0) 

{/it/} 



n 



(0) 



(3.17) 



(3.18) 



(3.19) 



denote the Lorentz-tensor symmetric and antisymmetric parts. These properties play an important role in ensuring 
transversality of the light-cone expansion, as will be shown in the following subsection. 



B. Transversality and twist— 3 operators 

Transversality of the twist-2 contribution. We now turn to the question how to maintain transversality (electro- 
magnetic gauge invariance) in the light-cone expansion of the current product. For the operator product, Eq.( |3.1| ), 
regarded as a function of the two points, x and y, current conservation implies that 

d d 
— = 0, —U^(x,y) = 0. (3.20) 



In terms of the center and relative coordinates, Eq.(3.2), these conditions take the form 



U { ^ } (z\X) = -—- n [H (*|X), (3.21) 



dzp 1 ' 2 8X^ 



d n [tlv] {z\X) = \-£-n {tlv} {z\X). (3.22) 



dz„ L ^ 1V 1 ' 2dX„ 



One sees that the transversality conditions relate the Lorentz-tensor symmetric and antisymmetric parts of the 
operator product; cf. Eq.( [2.8| ) for the Compton amplitude. 

In DIS the twist-2 part of the light-cone expansion, Eq.( |3.16] ), gives the leading contribution to the structure 
functions at large q 2 . In this case one is dealing with forward matrix elements of the twist-2 operators, so that all 



total derivatives of operators (i.e., derivatives d/dX) have zero matrix elements. The two conditions, Eqs.(3.21) and 
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( |3.22 ), decouple, and the properties Eq.(3.18) are sufficient to guarantee transversality of both the symmetric and 
antisymmetric parts separately. However, one immediately sees that in the case of non-forward matrix elements, 
as e.g. in DVCS, the twist-2 part alone is not transverse. This is because in the non-forward case total derivatives 
of operators generally have non-zero matrix element, and one can easily verify that in general 



JL n (0) 



* 0, 



JL n (0) 

dX u IH 



± 0. 



(3.23) 



As will be shown in Section IV, hadronic matrix elements of these total derivatives for DVCS kinematics are 
generally of order unity, i.e., not suppressed by factors m 2 j q 2 or t/q 2 . Thus, the violation of transversality by the 
twist-2 contribution alone cannot be neglected at the accuracy we are aiming for. Since the total operator product 
obviously is transverse, our conclusion can only be that in order to maintain transversality to the required accuracy 
we have to include operators of twist > 2 in the light-cone expansion. 

Including higher-twist operators. It turns out that the higher-twist operators needed to maintain transversality 
are contained already in the simple "handbag" term of the string operator expansion, Eq.(|3~5|); the subleading terms 
in 1/z 2 in the string operator expansion are not needed for this. The necessary operators are total derivatives of 
twist-2 operators, which are algebraically of twist 3. They naturally appear when taking into account the "rest" of 
the original string operator, which was dropped in taking the twist-2 part, and making use of the QCD equations 
of motion. 




In order to account for the twist > 2 parts of the string operators in Eq.(3.5) we first note that the differences 
between the symmetrized operators, E q.(|3.14| ), and the full string operators can be expressed in the form (for 
brevity we put X = in the following) [fjT 



^(-z/2) 7 ^(z/2) - [^(-z/2) 7 ^(z/2)] 6 



= / dt z a 
Jo 



—i>(-tz/2)~/ a iP(tz/2) - (a 



and similarly for the operator with Dirac matrix 70-75- To show this one uses that for any function, f(tz). 



d d 
z a — f(tz) = t—f(tz). 



dz a 



dt J 



(3.24) 



(3.25) 



The derivatives of the string operator in Eq.( 3.24 ) can be computed with the help of Eq^EiJ^). Of the terms with 



covariant derivatives of the fcrmion fields at the end points one vanishes by virtue of the QCD equations of motion, 



■ ■ ■ la Va 1p(x) = 0, 



= 0, 



(3.26) 



the other, with the help of Eq.(3.9), can be traded for a derivative with respect to the center coordinate, plus 
additional quark-gluon string operators. In this way one arrives at |13| 



^(-2/2)7^/2) - [^(-z/2) 7ff V(«/2)] 
r 1 



sym 



2 €aa ^ Za dX pJa 



dtt ip{-tz/2)j^ 5 ip(tz/2) 



string operators of type ip . . . F . . . tp. 



(3.27) 



A corresponding relation holds for the operator with Dirac matrix 7^75 on the L.H.S. and the operator with 7 7 on 
the R.H.S.: 



VK-z/2) 7( x75^/2) - ^(-2/2)7,75^(2/2)] 
i d f 1 - 

+ string operators of type %j) ...F .. .if). 



sym 



(3.28) 
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The precise form of the quark-gluon operators appearing here need not concern us; it turns out that these parts 
do not play a role in restoring transversality of the light-cone expansion. Note that Eqs.(3.27) and (3.28) relate 



the vector and axial vector string operators, which is possible with the help of the totally antisymmetric e-tensor. 



We emphasize that the operators appearing under the total derivative on the R.H.S. of Eqs.(3.27) and (3.28 ) are 
still full string operators, which have as yet no definite twist. In particular, they may again be decomposed into a 
symmetric (i.e., twist— 2) part and total derivatives, and so on. In fact, Eqs.(3.27) and (3.28) can be regarded as 



recurrence relations, repeated use of which allows to express the original string operator as the sum of its symmetric 
part and an infinite series of total derivatives of symmetric operators of arbitrary order (and quark-gluon operators, 
which we do not consider explicitly). We shall show below that these recurrence relations can be solved, and that 
the series of total derivative operators can be summed up in closed form. 

First-order correction in total derivatives. In order to understand how the total derivatives appearing in 
Eqs.([[27]) and ( |3~28| ) come into play in restoring the transversality of the light-cone expansion, it is instruc- 
tive to first take a look at the first-order correction in total derivatives to the twist-2 part of the operator product. 
To first order in d/dX we have (we do not write the contributions from quark-gluon string operators any longer) 



<K-z/2) 7 ^(z/2) = [^(-z/2) 7ff V(*/2)] 

d 



+ 2 e ' 7a ^ Za dXf 3 Jn 
+ terms of order 



dtt [ 

d d 
dXdX : 



-tz/2) lllb ij{tz/2)Y 



(3.29) 



and similarly for the operators with j a — * 7^75, 7^75 — > 7-,. When substituted in Eq.(3.5), the first-order correc- 
tions to the twist-2 string operators produce a contribution to the operator product which we denote by IL^J : 



n 



(1) _ z p J ■/ , \ 9 

2 4 1 l \ z ti9<?v 1 z uga^)£apP-i av 
TT Z OA 



fit/ 



Jo 



dtt [^)(-te/2)7 7 7 5 V'(te/2) + (z^-z)] s 



dx 



Jo 



dtt [ip(-tz/2)^ilj(tz/2) - (z^-z)Y 



(3.30) 



Let us see what the inclusion of this first-order correction in d/dX to the current product implies for the transver- 
sality conditions, Eqs.(3.21) and (3.22). By direct calculation one can easily show that 



1 d n ( 0) 



An« =I^_ n (0) 

dz„ [""I 2dX„ 



(3.31) 
(3.32) 



The equalities here are understood up to terms involving contracted total derivatives acting on the symmetric 
string operator, 



_d d_ 

dX^dX^ 



$(-z/2) 7 „iP(z/2)] 



sym 



or "total divergences" of the symmetric string operator 

d 



3X 



$(-*/2)7«^C*/2)] 



sym 



(3.33) 



(3.34) 



The hadronic matrix elements of derivatives of type Eq. ( p. 33] ) will be proportional to the square of the momentum 
transfer, t, so their contribution to the non- forward Compton amplitude will be suppressed by t/q 2 . Similarly, 



derivatives of type Eq.(3.34) will have matrix elements proportional to either m or t, and thus again lead to 
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power-suppressed contributions to the Compton amplitude. P| 
To verify Eq.(^3l|) one first shows that 



dx 



dz 



< I dtt [^(-te/2)7 7 7 5 V>(*z/ 2 ) 
'o 



l sym 



(3.35) 



and notes that the derivative z a (d/dz a ) acting on the string operator of argument tz can be replaced by t(d/dt), 
cf. Eq.(3.25). The R.H.S. of Eq.(3.31) is then obtained after integration by parts in the parameter, t. To verify 
Eq.(3.32) one has to make use of the fact the fact that, because of Eq.([3.14|), the symmetric part of the string 



operator satisfies 



jL[^-tz/2) l7l ^(tz/2)] sym = l-[^tz/2) ltll ^{tz/2)\ 



dz~, 



(3.36) 



and drop contributions of the type of Eqs.(3.33) and (3.34) 



Restoring transversality order-by-order in total derivatives. Is the inclusion of the first-order correction in total 
derivatives, Eq.( |3.30 ), in addition to the twist-2 contribution sufficient to restore transversality? Substituting 



n (°) 



(3.37) 



into the transversality conditions, Eqs.(3.21) and (3.22), and making use of the results Eq.(3.18) for the twist-2 



part and Eqs.(3.31), Q3.32 ) for the first-order correction, we find 



d 
dz^ 

d 
dz^ 



n 



n 



(0) 


{flu} 


(0) , 
[/*"] 4 


-np; 



1 d 



2dX IJ 
1 d 



n 



(0) 



2 dX b 



n 



[/XI/] 

(0) 



n 



(i) 



n 



(1) 



2dX„, M 



1 d 



n 



(i) 



'2dX u ^ 



(3.38) 
(3.39) 



Although the R.H.S. is of second order in total derivatives, it is not suppressed in the approximation we are 
considering, i.e., it is not proportional to derivatives of type Eqs.(3.32) or (3.34), which give power-suppressed 
contributions to the Compton amplitude. In fact, one can easily see that the "additional" total derivative in the 
operators on the R.H.S. appears in the combination z^d/dX^), which under the hadronic matrix element becomes 
—i(rz) and is not suppressed in DVCS kinematics. Thus, compared to the twist-2 part alone we have merely 
shifted the violation of transversality to terms one order higher in z^(d/dX^) = —i(rz). We conclude that, in 
order to obtain a transverse result, we have to include the total derivative operators to all orders: 



n 



transverse TT^ 



{IV 



n 



(i) 

//// 



n 



(2) 
Hi> 



(3.40) 



where the superscript (n) denotes the order in total derivatives. These contributions are generated by successive 
use of the recurrence relations, Eqs.( [3.27 ) and (3.28), which allow to express the "rest" of a symmetrized string 
operator in terms of total derivatives and quark-gluon string operators. (The latter turn out not to play a role in 
restoring transversality, so we shall drop them.) The transversality conditions connect the term of order n in total 
derivatives with that of order n+ 1; i.e., instead of Eq.( |3.31 ) and ( [3.32 ) one now has 



"For m 2 — and t = the matrix element of the twist-2 vector string operator has only contributions proportional to 
the momenta, p a and r a . Terms proportional to z a , which originate from z 2 -terms in the contracted operator, Eq.(3.14), 
come with factors of m 2 or t. 

^We stress again that the "order" in total derivatives here does not imply suppression of the contributions in t/q 2 . 
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9 n (n+l) = 1 9 



-n 



(n) 



2 dXf, ^ 



d n (»+i) 
3^ ^ 



1 9 rr(") 

2<9x M w 



(3.41) 
(3.42) 



All-order summation of total derivatives. Iterating Eqs.(3.27) and ( 3.28) ) , we express the string operators appear- 
ing in Eq.( ft.5| ) as the sum of its symmetric part, Eq.( |3.14| ), and an infinite series of total derivatives of symmetric 
string operators. This series can be summed up in closed form. The details of the calculation are presented in 
Appendix |A|. The result for the "fully deconstructed" string operators takes the form: 



V>(-2/2) 7ct V(z/2) 
= [^(-z/2) 7 ^(z/2)] sym 



+ 



8 



dX a 



i _d_ 



d 

dX 
i 

du u cos 



Jo 



du u sin 
d 

'8X 



iu 
~2 



dX 



[^(-«2/2)7 7 ^(u«/2)] 



sym 



-mz/2)7 7 y 5 i/;(uz/2)Y 



contracted operators of the type (3.33) or ( ft .34 ) 
quark-gluon operators, 



(3.43) 



and similarly for the operators with ^ a — » 70-75, 7« 75 — > 7«- [Note that when expanding to first order in d/dX we 
recover Eq.(3.2£).] We shall refer to the part of the string operator involving only symmetrized operators and their 
total derivatives as 



[^(-z/2) 7 ^(z/2)] 



sym + total der. 



(3.44) 



In Eq.(3.43) we can substitute on the R.H.S. the explicit representation for the symmetric string operators, 
Eq.( |3.14 ). Performing an integration by parts over u in the second term of Eq.( 3.43 ), and using the identity 
(here v = ut) 



du I dt f(u) g(ut) 



du / dvf(u)g(v) 



dvg{v) 



duf{u) 



(3.45) 



to simplify the third term, we obtain 

^(-z/2)7^/2)] sym + totaldcr - 



= I dv < cos 



IV 

i d d 



d 

'OX 



d iv 
— + —sin 

oz a 2 



dx 



dX a 



ip{-vz/2)zip(vz/2) 





iu 




/ dv du cos 




(4)] 


/o Jv 


~2 





ip(—vz/2)zj5ip(vz/2). 



(3.46) 



The parameter integral in the second term can be carried out and gives 

d 



du cos 



1 u 
~2 



dx 





iv 




sin 




(4)] 


2 





d 

'dX 



(3.47) 



An analogous formula applies to the operators with j a — > 70-75; on the R.H.S. one should replace I — > £75, £75 — > z. 
Note that to zeroth order in total derivatives, the R.H.S. of Eq.( 3.46 ) just reduces to the usual symmetrized operator, 
Eq.(p^). 

Finite translation of symmetric string operators. The sine and cosine functions involving the total derivative, 
appearing in Eq.(3.43), are nothing but the odd and even part of the finite translation operator 
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exp 



±V 2{ Z dX) 



(3.48) 



which affects a translation of the center of the string operators by a distance ±vz/2. Taking this into account, we 
may alternatively write the result for the string operator containing total derivatives of symmetrized operators, 
Eq.( [3.46| ), in the form 



-z/2) 7ct ^(z/2)] 



sym + total der. 



dv 



d 
dz a 

d 



v d 



2dX 
v d 



ip(-z/2)z i>[{2v - 1)2/2] 
$[-(2v- l)z/2] zip{z/2) 



4 OX/3 OZ-y j 



dz a 2 dX a / 

d d f 1 f 1 

dv / du iip[—(v + u)z/2] zj5ip[(v — u)z/2] 



+ ip[-(v - u)z/2] zj 5 ip[(v + u)z/2] 



(3.49) 



We see that taking into account total derivatives of symmetrized operators to all orders in the decomposition of the 
string operator amounts to certain finite translations of the symmetrized operators, along the direction defined by 
the separation of the fields in the original operator, z. It is interesting to note that the endpoints of the symmetrized 
string operators under the parameter integrals lie inside the interval [— z/2, z/2] defined by the original operator; 
this can be regarded as a consequence of locality. 

Transverse extension of leading light-cone singularity. Substituting the closed-form expression for the decompo- 
sition of the string operator in symmetric operators and total derivatives thereof into Eq. ( |3.5| ) we can now obtain 
the expression for the transverse extension of the twist-2 part of the light cone expansion, Eq.(3.40): 



n 



transverse 



4gj { s, pva $(-z/2M(z/2) - (z -> -z)] sym + total d ° r 

r 7/ /n\ 1 i ic\\ / \1 s Y m + total der 

-te^pua [tp(-z/2)^ (T ^ 5 2p(z/2) + (z — > -z)\ " 



} 



(3.50) 



One can verify that this expression satisfies the transversality conditions, Eqs.( 3.21 ) and ( 3.22 ), up to terms 
involving operators of the type Eqs.( |3.33| ) or (3.34), which give power-suppressed contributions to the Compton 
amplitude. [To see this, one needs to make use of the identity Eq.( 3.25 ) and integrate by parts over the parameter 



v.] We do not bother to make the string operators traceless here. As explained above, z -terms in the operators 
need to be included explicitly only when keeping corrections of order t/q 2 and m 2 /q 2 to the DVCS amplitude. Note 
that in this case one would need to include also operators of the form Eqs.( 3.21 ) and ( |3.22] ) in the decomposition 
of the string operators, Eq. (|3.43 ). 

We can substitute in Eq.( 3.5C ) the explicit expressions for the vector string operators in terms of the contracted 



string operators, Eq.(3.46), and express our result for the transverse generalization of the leading term light-cone 
expansion directly in terms of the contracted operators. Grouping together contributions from the operators with 
Dirac structures z and £75 we obtain 



Trtransvcrsc 



7T Z Jq 



til- ^ COS 
1 

2 



d 



iv ( d 
~2 VdX 



d 



du cos 



iu ( d 
~2 VOX 



_d d_ 

dXp dz 1 



^We are grateful to V.M. Braun for bringing this point to our attention. 
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x [%p(-vz/2)z^(vz/2) + (z -> -z)] 



dv 



-i cos 

i 
2 



Y Vdx 



_d_ 

' dz a 



v . 
— sin 
2 



iv 




2 





cos 



y rax 



_d 9_ 



x vz/2)z"f5i^(vz/2) — (z — » — z)l 



(3.51) 



This form is useful because the two parts involving the operators ^>z^> and i/jz'ys'ip turn out to be individually 
transverse to the accuracy stated above. (Note that both these terms contain a symmetric and an antisymmetric 
tensor part.) This is necessarily so: There is in general no dynamical relation between the matrix elements of 
the operator ipzip and those of the operator ipz'ysip (in the forward limit they are related, respectively, to the 
unpolarized and polarized quark distributions in the hadron), so their contributions to the light-cone expansion 
must be individually transverse. 



C. Quark matrix elements: Recovering the handbag diagram 

Before turning to the hadronic DVCS amplitude, it is instructive to compute the matrix element of the leading 
term in the light-cone expansion, Eq.( 3.50| ), between quark states. This simple exercise serves two purposes. First, 
we can verify that the Compton amplitude for a quark target, obtained from Eq.(3.50), is — up to terms of order 
t/q 2 — equivalent to the result of the simple "handbag" graph in a theory of free quarks. This may not come as 



a surprise, as the operator result, Eq.(3.5C), was obtained by way of a reduction of the simplest two-point string 
operator, whose graphical analogue in the free theory would just be the Compton amplitude for the free quark. 
Still, it is worthwhile to check explicitly that the operators of twist > 3, which were dropped in the process of 
reduction of the string operator, are not neeeded to reproduce the "handbag" result. Second, when computing 



the quark Compton amplitude from Eq.(3.50) we shall see that characteristic singularities appear in terms with 
denominators quartic in momenta, which are a consequence of the finite translations of the string operators in the 



gauge-invariant expression, cf. Eqs.(3.46) and (3.49). In the free quark case these singularities cancel when adding 
the contributions from the two twist-2 operators, ipzip and ipz^ip. Similar singularities we shall also encounter 
in the DVCS amplitude for a hadronic target; however, in this case they do not cancel completely, leaving certain 



divergences in the tensor amplitude for transverse photon polarization (see Section IV). 

Quark matrix elements of symmetrized string operators. Let us compute the matrix element of Eq.( [3~5C ) between 
free massless quark states. Their wave functions are plane waves of momenta p ± r/2, multiplied by Dirac spinors 
satisfying 



I U = 0, 



U (£--] = 0. 



(3.52) 



The quark matrix elements of the totally symmetric (twist-2) light-ray operators (with center coordinate X = 0) 



arc 1 



{p-r/2\i>(-z/2)\ Z U(z/2)|p + r/2) = U I \ue~^ z \ 



(3.53) 



Only the average momentum, p, enters in the phase factor containing the separation of the fields, z, in accordance 
with translational invariance. From these one can easily construct the matrix elements of the vector and axial 



vector string operators including symmetrized operators and their total derivatives, Eq.(3.46). We find 
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(p - r/2| ^(-z/2) 7 ^(z/2)] 



sym + total der. 



\P + r/2) 



+ / dve~ m{pz) \ ivcos[v{rz)/2\U[{pz) lrJ -p a z\U 



■ s'm[v(rz)/2] U [(rz)j a — r a z] U 



+ ^aaPjZar^ j du cos[u(rz)/2] U (7 7 - ivp 1 z)^ 5 U 



(3.54) 



Here the first three terms originate from the contracted operator ijjzip, the fourth from ipzj^ip. We have performed 
an integration by parts in the integral over the parameter v in the term proportional to Uj a U, in order to isolate 
the first term, which coincides with the free quark matrix element of the vector string operator, iJ;(—z/2)'y ( jijj(z/2). 
Note that all terms but this one vanish in the limit of zero momentum transfer, r — > 0. In order to simplify the 
second and third term we make use of the following identities for quark bilinears, which follow from the Dirac 
equation for the quark spinors, Eq.( 3.52] ), and the identities for products of three gamma matrices: 



- U [(rz)j a - r„z\ U = i<L a<x pi z <xV& U^ 5 U 



(3.55) 
(3.56) 

(similar identities hold with 7 replaced by 775 everywhere). Furthermore, to simplify the third term, we perform 
an integration by parts in the piece with 7 7 and make repeated use of Eqs.(3.55) and (3.56): 

i 

= -^aap^s^ZaZgrpr,: U^qU 
% — 

= (rz) U [(rz)7„ — r v z\ U + terms oc r 2 



-(rz)e aa f3 1 z a pi3 C/7 7 7 5 t/. 



(3.57) 



Using these identities one can easily convince oneself that the three terms in the square bracket in Eq.( |3.54| ) add 
up to zero. The result for the vector string operator including symmetrized operators and total derivatives is thus, 
up to terms of order r 2 — t which were dropped in Eq.(3.57), completely given by: 



(p-r/2\ [^(-z/2) 7 ^(z/2)] 



sym + total der. 



\p + r/2) = e -^) U la U, 



(3.58) 



which coincides with the matrix element of the vector string operator between free quark states. A similar result 
is found for the axial vector string operator: 



(p - r/2\ [<H-z/2) 7CT 75^/2)] Sym + ^ ^ \ P + r/2) - U lalb U. 



(3.59) 



Compton amplitude for the free quark. In particular, Eqs.( 3~58| ) and (3.59) imply that the Compton amplitude 
for a free quark, calculated from the transverse generalization of the leading term in the light-cone expansion, 
Eq.(3.5C), is up to terms of order t/q 2 identical to the "handbag" Compton amplitude: 



d 4 ze l{qz) (p-r /2\Il t ™ 



transverse 



\p + r/2) = U 



(p — q) 2 + iO 



+ (<? -» -Q,IJ>*-> v) 



U. 



(3.60) 



This amplitude is known to be transverse, by virtue of Eqs.(3.52). To summarize, our calculation shows that the 
operators dropped in Eq.( |3.50]) indeed contribute to the free quark Compton amplitude only at power-suppressed 
level. 
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It is instructive to take a look at the contributions of the individual terms under the parameter integral in 
Eq.( [3.54 ) to the Compton amplitude. In the first and second term in the bracket, the sine and cosine functions of 
argument v(rz)/2 can be combined with the overall exponential factor to give exponentials of the form 



e -i[v(pz)^fv(rz)/2] 

The contribution of these terms to the virtual Compton amplitude then leads to integrals of the form 



dv v Id ze 



i(qz) 



2z ^ p QL 



-i[v(pz)^v(rz)/2] 



dv 



[{q - vp ± vr/2) 2 + iOp 



_7r 2 (z 2 -i0) 2 

In the generalized Bjorken limit these denominators become (up to terms of order t/q 2 ) 

(q — vp± vr/2) 2 — 2(pq)(—v ± vr\ — £), 
where £ and r\ are defined in Eq.( |2.9| ). In the DVCS limit 77 = £, and thus 

(q- vp±vr/2) 2 = 2{pq) { 

In the case of the upper sign one is led to integrals of the form 

dv 
v + iO' 



x tensor. 



(3.61) 



(3.62) 



(3.63) 



(3.64) 



(3.65) 



which diverge logarithmically. This singularity is canceled by a corresponding singularity produced by the third 
term in the bracket in Eq.( 3.54 ), which represents the contribution of the twist-2 operator with Dirac matrix £75, 
as we have shown above, using the equations of motion for the quark spinors. As we shall see in Section IV C, in 
the case of hadronic matrix elements, where the quarks are not on mass shell, this cancellation takes place only 
imperfectly, and a real divergence occurs in the Compton amplitude for transverse photon polarization. 



IV. DVCS AMPLITUDE INCLUDING KINEMATIC AL TWIST-3 TERMS 
A. Spectral representation of twist— 2 matrix elements 

We now use our result for the transverse generalization of the leading term in light-cone expansion of the product 
of electromagnetic currents to compute the DVCS amplitude for a hadronic target. In this way we shall obtain an 
approximation for the DVCS amplitude which is manifestly transverse up to terms of order t/q or m 2 /q 2 . 

In order to compute the hadronic Compton amplitude we need to supply parametrizations of the hadronic matrix 



elements of the basic symmetrized (twist-2) string operators in Eq.( 3.50 ). For simplicity we restrict ourselves to the 
case of a pion target, which has spin zero and simple charge conjugation properties. Specifically, we shall consider 
the isoscalar component of the DVCS amplitude for the pion, and suppress the overall factor resulting from the 
quark charges. We shall neglect the pion mass and put m = 0. The generalization of the expressions presented 
below to higher-spin targets or to other isospin components of the amplitude is straightforward. 



z 2 -dependence. When computing the Compton amplitude, Eq^^l), from the light-cone expansion in coordinate 
space, the integral over the separation of the fields, z, runs over the whole 4-dimcnsional space. Thus, in principle 
we need to provide a parametrization of the matrix elements valid everywhere in z, not only on the light cone. 



Consider the hadronic matrix element of the contracted string operator, Eq.(3.14), for a massless spin-0 hadron 



.,2. 



On general grounds it can be regarded as a function of the variables pz, rz and tz 

(p-r/2\ip(-z/2)zip(z/2)\p + r/2) = function (pz, rz; tz 2 ). (4.1) 
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The dependence on z 2 is through the combination tz 2 ; the reason being that t is the only dimensionful invariant in 
the problem. Since we shall drop terms of order t/q 2 in the Compton amplitude, it follows that we can effectively 
neglect the z 2 -dependence of the matrix element in the integral over z. In this sense it is sufficient to provide a 
paramctrization of the matrix elements of the operators only on the light cone, z 2 = 0. 

One would have to take into account the z 2 -dependence of the matrix elements if one wanted to compute 
kinematical corrections to the DVCS amplitude of order t/q 2 or m 2 /q 2 . Such corrections are analogous to the 
well-known target mass corrections to DIS [^9|,[30| . In this case one would have to remember that the definition of 
the twist-2 part of the contracted operator implies that it satisfies the harmonic condition [cf. Eq.( 3.15| )1 



□,[^(-z/2)^(z/2)] twist " 2 = 0, (4.2) 



which is equivalent to the condition that the local operators in the Taylor expansion in z be traceless. Eq.(4.2) 
allows to reconstruct the z 2 -dependence of the matrix element of the operator from "initial conditions" given on 
the light-cone, z 2 = 0. In our calculation here we shall explicitly put t and m 2 to zero, so we do not need to worry 
about trace subtraction. In this sense we do not need to distinguish between the symmetrized string operator 
Tp(—z/2)zip(z/2) and its twist-2 part. 

Spectral representation for the matrix elements. The most convenient way to parametrize the hadronic matrix 
elements of the symmetrized string operators is by way of a spectral representation, i.e., a decomposition in plane 
waves. In the case of the pion the matrix element of the pseudoscalar operator, ip(— z/2)z , y5ip(z/2), is identically 
zero. (The forward matrix element of this operator would be parametrized by the polarized parton density.) Thus 
we need only the parametrization of the matrix element of the contracted scalar twist-2 operator. For the isoscalar 
component of the matrix element it is of the form 

(p - r/2$(-z/2)zi/j(z/2)\p + r/2) 

= 2{pz) / dx / dae- l(kz) f(x,a) 

J-l J-l+\3i\ 

+ (rz) J da e~ la{rz)/2 D{a) 

+ z 2 -dependent terms, (4.3) 

where 

k = xp + —r. (4.4) 

In the partonic language, k is the average of the outgoing, xp + (1 + a)r/2, and incoming, xP — (1 — a)r/2, parton 
momenta. Here, f(x,a) is the so-called double distribution |p|, pT| , p2[ , while D(a) is the Polyakov- Weiss (PW) 
distribution amplitude absorbing the (pz)-independent terms of the matrix element ||. They obey the symmetry 
relations 

f(x,a) = -/(-£, a), (4.5) 
f(x,a)= f{x,-a), (4.6) 
D(a) = -D(-a). (4.7) 



Eq.(|4.q) is the general symmetry of double distributions noted in Ref. pi] , while Eq.(4.£) is specific to the their 
C-even components. The z 2 -dependent terms in Eq.(4.3) are in principle calculable in terms of f(x,a) and D(a) 



by solving the harmonic condition, Eq.( |4.2| ), to all orders in z 2 . However, since each z 2 is accompanied by a factor 
of either t or m 2 , we can neglect them, as explained above. 
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B. Matrix elements of vector— type string operators 



From the parametrization of the matrix element of the contracted string operator, Eq.( |4.3| ), we can obtain the 
matrix elements of the the vector and axial vector string operators, Eq.( 3.46 ), including the kinematical twist-3 
contributions represented by the finite translations. Let us consider first the part of the matrix elements coming from 
the double distribution term in Eq.(4.3); the contributions from the PW-term can be included later. Substituting 
Eq. (fL~3|) into Eq.( 3.46 ), and using the fact that under the matrix element the total derivative of string operators 
turns into the momentum transfer, 



(p — r/2\i 



d 



\p + r/2) = r a (p-r/2\...\p + r/2), 



(4.8) 



we obtain: 



(p-r/2\ [iP(-z/2) la iP(z/2)] 



sym + total der. 



\p + r/2) 



dx 



l-\x\ 



—H-|3 



daf(x,a) / dve~ m{k ^ 



x {2 [p a — ivk a (pz)] cos[v(rz)/2] + vr a (pz) sin[v(rz)/2]} 



dx 



l-\x\ 



1+1*1 



daf(x,a) [2p a e-^ 



dvve lv ( kz > {sin[w(rz)/2] [p a (rz) - r a (pz)} - 2i cos[v(rz)/2] \p a (kz) - k a (pz)]} 



(4.9) 



In the last step we have performed an integration by parts in the integral over the parameter v in the term 
proportional to p& • As described in Section ||, in DVCS kinematics we can express the momentum transfer in the 
form [cf. Eq.|Tl5|)] 



2£P + A, 



(4-10) 



where A is of order |i| 1//2 . In the partonic language the term proportional top would be the longitudinal component, 
and A the transverse component. Similarly, we write the "active" quark momentum in the form k = (x + £,a)p + 
aA/2. Noting that in the square brackets in Eq.(4.9) only the A-parts of the momenta r and k contribute, 



p a (rz) - r a (pz) = p a (Az) - A a (pz), 
p a {kz) - k a {pz) = - [p CT (Az) - A a {pz)\ 



(4.11) 



we obtain for the matrix clement 



(p-r/2 | [i,{-z/2)^{z/2)] 



sym + total der. 



72) 



d.i: 



1+1*1 



da f(x, a) 2p a e 



,-i(x+£a)(pz)-ia(Az)/2 



+ [A a (pz) -p a (Az)]J^ dvve -Mx+ia)( pz )-w a (A z )/2 { sin [ 5 ( rz )/ 2 ] -iacos[v{rz)/2}}^ . 



(4.12) 



Introducing skewed distributions. Since we aim to compute the DVCS amplitude up to terms of order t/q 2 , we 
can expand the matrix elements of the vector string operator in A ~ l^ 1 ^ 2 - This will allow us to express the matrix 
elements in terms of "skewed" distributions, i.e., one-variable spectral functions depending on the kinematical 
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variable £ as an external parameter. Expanding the exponential factors and the sine and cosine functions in 
Eq.(p^) in (Az) we obtain, to first order in Ar 



<p - r/2\ $(-z/2M(z/2)] 



sym + total der. 



dx I da f(x, a) ( 2p a 

1 J-l + \x\ 



1C( , . 

1 (Az 

2 V 



\p + r/2) 

—i{x+£a)(pz) 



(4.13) 



+ [A a {pz) - Pa (Az)}J^ dv v e~ m{S:+ t a){pz) {sm[v£(pz)]-iacos[vZ{pz)}}) . 

Now the spectral parameters, x and a, appear in the exponential factors only in the combination 

x = x + £a, 

so the information entering in the matrix element is effectively contained in certain one-dimensional "reductions'' 
of the double distribution, fix, a). We introduce two skewed parton distributions: 



H(x,0 
A(x,0 

This is visualized in Fig. |l|. 



i ,i-|x| r i 

dx / da 5(x — x — £cv) f(x, a) < 
l J-i+\x\ I a 



(4.14) 





a 








1 










\ 1 i 

\ \ 

\ * 

\ 1 \ 




-1 




xy 


1 X 




-1 







FIG. 1. Graphical representation of the reduction formulas, Eq.(4.14), defining the SPD's, H(x,£) and A(x,^), in terms 
of the double distribution, f{x,a). The range of the arguments x and a of the double distribution is the "rhombus" defined 
by | J | + \a\ < 1. The integration is over the line defined by the condition x + £a — x and lying inside the range of the 
variables x , a; it is indicated by the fat line. The variable x defines the intersection of this line with the i-axis, the variable 
£ its deviation from the vertical. For x — £ the line runs through the upper corner of the rhombus. 



*** Since we plan to drop terms of order t/q 2 in the Compton amplitude it seems sufficient to expand the matrix elements 
to first order in A, since A 2 -terms would already be proportional to t. There is, however, a subtle point concerning 
contributions to the Compton amplitude of tensor structure A M A„, which one needs to include in order to convert an 
approximately gauge-invariant structure in the tensor amplitude into an exactly gauge-invariant one. Although the A M A„ 
term is formally of order 0(A 2 ), it can be obtained only from the contribution corresponding to the leading light-cone 
singularity, z/(z 2 ) 2 , cf. the discussion at the end of Section [VC. 
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As a consequence of the symmetries of the double distribution, Eqs.( |4.5| ) and ( |4.q ), the functions H and A satisfy 
the symmetry relations 



H(x,0 = -H(-x,0, H(x,0 = H(x,-0, 

A(x,Q = A(-x,0, A(x,0 = -A(x,-0- 

Furthermore, because of the antisymmetry of the combination a/(i, a) with respect both to x and a we have 



(4.15) 



dxA(x,£) = 0. 



(4.16) 



Note that the integral over x is taken over the < x < 1 interval. Hence, the distribution A(x, £) cannot be a 
positive-definite function on < x < 1. 

Combining in Eq.(4.13) similar terms we obtain the following representation of the matrix element in terms of 
the skewed parton distributions: 

(p-r/2\ [^(-z/2) 7 ^(z/2)] sym + totaldcr - b + r/2) 



dx ( 2p a e 



-ix(pz) 



H(x,0-^A(x,0 



+ \A a {pz) - Pa(Az)] j dvv e ~ lvx (P z > {H(x, f ) sin[uf (jw)] - iA(x, f) cos[v£(pz)]} \ . (4.17) 

The cosine and sine functions can be represented as the sum or difference of exponentials. Combining them with 
the overall exponential factor, e ~ lvxl J> z ) ; ne gets vx ± vt; combinations. However, using the symmetries in x of the 
SPD's, Eq.( 4.15 ), one can easily arrange that only vx + vt; appears. This finally gives 

<p - r/2| ^(-z/2) 7 ^(z/2)] sym + total dcr - \ P + r/2) 



dx 2p a e 



-ix(pz) 



H(x,0~^A(x,0 



- [A a (pz) - P(7 (Az)] I dvv 



i(vx+v£)(pz) , e i(vx+v£)(pz) 



[H{ X) Z)-A(x,0] 



(4.18) 



The axial vector string operator. In a similar fashion we compute the matrix element of the axial vector string 
operator, Eq.( 3.46 ). Substituting the double distribution part of the spectral representation of the contracted 
operator, Eq.((D|), into Eq.(|3~46|), we obtain 



(p-r/2\ [^(-2/2)7*75^/2)] 



sym + total der. 



\p + r/2) 



= 2e aal3l z a r0 / dx 



da f(x, a) [ dv [p 7 - ivk^pz)] e~ m ^ S ^MI3 



1+ X 



(rz) 



(4.19) 



Again integrating by parts in the p 1 term, and keeping only terms linear in the transverse component of the 



momentum transfer, A, this turns into 



(p-r/2\ $(-*/2)7a7feV'(«/2)] 



sym + total der. 



b + r/2) 



1-1*1 



1+1*1 



= Cvapy z a Ap Pl I dx I daf(x,a) / dvve ™( x +t a )(P z ) { C o$[v£(pz)] + ia sin[v£(pz)]} 



(4.20) 



One sees that this matrix element can be expressed in terms of the same skewed distributions, H and A, as the 
one of the vector operator: 

( P -r/2|[^(-z/2) 7(775 ^(z/2)] sym + totaldor - 



\p + r/2) 



\e aa(il z a A fi p 1 I dx I dvv 



e -i(vx+v£)(pz) _ e i(vx+v£)(pz) 



[H(x,0-A(x,0] 



(4.21) 
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Hence, the matrix elements of both the vector and axial vector string operators, including the twist-2 and 
kinematical twist-3 contributions, are completely described by the two skewed distributions, H and A, which, in 
turn, are determined by the original double distribution, Eq.( |4.3| ), through the reduction formulas, Eq.( 4.14 ). 

Contribution from the PW-term. In order to obtain the full matrix elements of the vector and axial vector 



operator we must add to Eqs.( 4.18 ) and ( 4.21 ) the parts coming from the PW-term in the spectral representation 
of the matrix element of the contracted operator, Eq.(4.3). Substituting this term in Eq.( [3.46| ), and going through 
the same steps as described above for the double-distribution part, one easily sees that the contribution of the PW 
term to the matrix element of the vector operator is of the form 



(p-r/2\ [<K- 2 /2) 7ct V(2/2)] 



sym + total der. 



\p + r/2)p W -Terro = 



dae- ia(rz)/2 D(a). 



(4.22) 



Note that the PW-term has a simple structure, corresponding to a parton picture in which the partons carry 
fractions (l±a)/2 of the momentum transfer, r. Since only one momentum, r, is involved, this term can contribute 
only to the totally symmetric part of the vector string operator, and thus "decouples" in the reduction relations, 



Eqs.( p.28| ) and (3.28). For the same reason the PW term also does not contribute to the matrix element of the 
axial vector string operator: The derivatives with respect to both coordinates, X and z, give rise to the momentum 
transfer, r, whence the contraction with the e-tensor gives zero. 



C. DVCS amplitude for pion target 

Having at hands the parametrizations of the matrix elements of the vector and axial vector string operators 
we can now compute the DVCS amplitude for a pion target. Substituting the parametrizations, Eqs.( 4~i~8| ) and 
(4.21), into Eq.( |3.5d ), and performing the Fourier integral over the distance, z, we obtain the hadronic Compton 
amplitude in the form 



= I dx 



d 4 z e i{qz) - 



2iz n 



X \ 2s [ipi/aPa 6 



-ix(pz) 



TT 2 ( Z 2 -iO) 2 



+ 2 Sppva [A a (pz) -p a (Az) 






x dvv 


e ~i{vx+v(,){pz) 


Jo 




i 






x dvv 




Jo 




(q -t-q,/**-* 


V). 



[H(x,S)-A{x,S)] 



iH(x,0-A(x,0} 



(4.23) 



Using the well-known identity for the contraction of two e-tensors one sees that the second and third term here 
can be combined to give two terms involving only exponential factors with positive and negative sign, respectively, 
which give identical contributions. The integral over z is readily computed, using the basic formulas: 

Z n In 



d 4 ze^ 



2n 2 (z 2 -i0) 2 l 2 + i0' 



d 4 ze 1 ^ ™? a 



11 n l 



2ir 2 {z 2 -i0y 



(i 2 + ioy 



(4.24) 
(4.25) 



Here I is to be replaced by the sum of all momenta appearing in the exponential factors, and I 2 in the denominators 
can be simplified using the special relations between invariants in DVCS kinematics, cf. Section [n]. Making use of 
the symmetry properties of the pion SPD's, H and A, Eq.(4.15), the result can be assembled in the form 
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Tpv = -^—^ | 2 [p^ + q^p v - g^, v {pq) + 2£p (tl p„] J dx 



x-^ + iO 



x-£ + iO 



J dx[H(x, S)-A(x,Z)] 



dv v 



- {vx + vQPft 
(£ + vx + v£ - iO) 2 



A, 



q v + (vx + v£)p v 



(4.26) 



We can bring this result into a form which is manifestly transverse with respect to both the incoming and outgoing 
photon momenta. By simple algebraic rearrangement, and making use of the integral identity 



dv 1 



-vx — v£ 



1 



(£ + vx + v£ - iO) 2 x + £-iO 



= dvi 



3£ 



(£ + vx + v£ - iO) 2 ' 



(4.27) 



which holds if x ^ ±£, we first write Eq.( 4.26 ) as a sum of three terms which are individually transverse JT(|. In 
addition, we then perform an integration by parts in x in order to remove the quadratic denominators. Here we 
use that 



H(x,0-A{x,0 



for 



±1, 



(4.28) 



which naturally follows from the reduction formulas, Eq.(4.14). In this way we finally obtain 
1 



T 



2(pg) 



PnQu + q^Pv - 9^u(pq) + 2£PnPv ~P^— + —Pv 



1 a g^iij 

dx 

_ x x-£ + iO 



(q, L + 3^pn)A„ / dx dv 



dv ■ 



R(x,0 



&n{<lv + &v) I dx j 
J -I Jo 



£ + vx + v£ — iO 



— £ + vx + v£ + iO 

where R(x, £) is a short-hand notation for the following combination of derivatives of the SPD's: 

R(x,0 = ^-H(x,0-^A(x,0- 



(4.29) 



(4.30) 



All three terms here are individually transverse up to terms of order t. The decomposition of the Compton 
amplitude here is identical to the one introduced in Ref. |l6|. The first term is the twist— 2 part with the tensor 
structure corrected as suggested by Guichon and Vanderhaeghen fl5|| . The second and third term represent twist-3 
contributions. Note that the second term contributes only to the helicity amplitude for a longitudinally polarized 
initial (virtual) photon, the third term for a transversely polarized one. 



rft 



One observes that the second SPD, A(x, £), enters in Eq.( 4.29 ) only through its first derivative, cf. Eq.( 4.3C ). 
From the reduction formulas, Eq.(4.14), one sees that the ^-derivative of the function A is related to the ^-derivative 
of H by 



' T Tn general the virtual Compton amplitude for a spin-0 target has five independent structures. The fact that in the 
approximation considered here (twist-2 plus kinematical twist-3, neglecting t/q 2 and m 2 /q 2 contributions) only three of 
these structures are non-zero is analogous to the absence of the longitudinal structure function in DIS at twist-2 level due 
to the Callan-Gross relation. 
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(4.31) 



so that instead of Eg. (4.30) one could equivalently write 



(4.32) 



One could thus express the amplitude, Eq.(4.29), completely in terms of the first SPD, H, and its ^-derivative. In 
this way one would obtain the form of the Compton amplitude given in Ref. |2(J . 

Inspecting the singular integrals. We now inspect the singular integrals which appear in the the expression for 
the Compton amplitude, Eq.( 4.29j ). The integrand in the first term has only a simple pole at x = £; this integral is 
of the same type as that obtained keeping only the twist-2 part of the matrix elements. The integral exists if the 
SPD, H, is continuous at x = £, which is the case for SPD's derived from DD's which are less singular than 1/x 2 



for x — > and otherwise continuous 



The second and third term in Eq.(4.29) involve additional parameter 
integrals, so their convergence needs to be investigated carefully. In the second term, which contributes to the 
hclicity amplitude for longitudinal polarization, we can perform the parameter integral over v, excluding the value 
x = £ for the moment: 



dv 



1 



£ + vx + vl; — iO 



hi 



x + £ - iO 
2£-i0 



(4.33) 



This function has a logarithmic singularity at x = — £, while it is regular at x = £. The integral over x in Eq.( 4.29] ) 
exists if R(x, £) is bounded at x = — £, which again is the case in the DD-based models described in Ref. Q. We 
thus obtain a finite answer also for the second term contributing to longitudinal photon polarization. 



Let us now turn to the third term of Eq.(4.29). Here we are dealing with an integral of the type (remember that 
V = 1 — v) 



dv 



—l;+VX + v(;+ iO 



dv 



o v(x-£) + iO 



(4.34) 



This integral exhibits a logarithmic divergence at the lower limit. One may hope to get a finite result only if the 
integral 



no = 



\x- R ^ 



£ + i0 



(4.35) 



vanished. However, one can easily convince oneself that this quantity is generally non-zero. In fact, representing 



i?(x,£) as in Eq.(4.32), and performing an integration by parts in the term with the rr-derivative, one finds that 



7 «> - 5 



iO 



(4.36) 



which is the ^-derivative of the twist-2 contribution to the DVCS amplitude. Since the latter generally has a 
nontrivial ^-dependence, determined by the shape of the SPDs, we conclude that the parameter integral Eq.( 4.34 ) 
really is divergent. Thus, in our approximation the twist-3 part of the tensor amplitude T M „ corresponding to 
transverse polarization of the initial photon contains a divergence. A similar observation has recently been made 
by Kivel et al. @. 

Absence of divergences in the physical DVCS amplitude. One may wonder if the divergence occurring in the tensor 
amplitude for transverse polarization of the initial photon affects the physical DVCS amplitude. Fortunately, this is 



not so. One can easily see that the tensor structure of the third term of Eq.(4.29) involving the divergent integral 



"''Continuous SPD's were also obtained in model calculations of the SPD's at a low scale in the instanton vacuum [B31 
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(4.37) 



is just a truncated (to first order in A) version of the exactly gauge invariant combination, 

A„q 2v = A^ + rv/2) = A^^ + ^ + A^/2), (4.38) 

which has zero projection on the polarization vector of the final real photon, ti v : (£212) = 0. In fact, the structure 
A )1 g2 I/ is what one would obtain if one calculated the DVCS amplitude directly from the original "unexpanded" 
form of the DD parametrization of the matrix element, Eq.( 4.12] ). It appears from the term with the exponential 
factor of the argument 

- i[v(x + £a)(pz) +v£(pz) + (va + v)(Az)/2], (4.39) 



which is obtained by combining the sine/cosine functions and the exponential factor in the second term of Eq.( 4.12] ) 



In the Compton amplitude, this term gives rise to a contribution in which the argument of the quark propagator is 

q + v{x + £a)p + v£p + {va + v)A/2. (4.40) 



Since (Aq), (Ap) and A 2 are negligible, the denominator factors in Eq.( 4.29 ) remain the same as above. In the 
numerators, however, writing (va + v)A/2 as [1 — (1 — a)v]A/2, we observe that the A fl (q L , + term of the above 
result converts into A M (q v + ^p u + A u /2) — A^q 2vl plus a contribution of the form vA^A u . The latter corresponds 
to a "new" SPD defined by a reduction formula of the type of Eq.(1.14) but with integrand (1 — a) 2 f(x, a). Due 



to the extra factor v, the u-integral for the latter contribution is finite. Thus, inclusion of the A^A^ term in the 
matrix elements and the Compton amplitude relegates the divergence in the tensor amplitude to the (exactly) 
unphysical part which does not contribute to the photon helicity amplitudes. 

It may appear that retaining the A M A„ term in the Compton amplitude would exceed the accuracy of our 
approximation, since this term is formally of order t, and thus of the same order as corrections from z 2 -terms in 
the matrix elements which we did not take into account. However, the A^A^ term is distinct in the sense that it 
can only arise in the way described above, i.e., from keeping finite A in the exponents of the matrix elements, while 
z -terms in the matrix elements can give rise only to corrections of order t/q 2 to already existing tensor structures. 
In this sense the inclusion of z 2 -terms would not spoil the above argument proving that the physical Compton 
amplitude is divergence-free up to terms of order t/q 2 . 

To summarize, for the physical DVCS amplitude (i.e., the photon helicity amplitudes) we find no evidence for 
factorization breaking in the kinematical twist-3 contributions, both in the \Jtjq 2 and t/q 2 terms. It is quite 
possible that factorization breaks down at the t/q 2 level, but in order to establish this one would need to analyze 
the contributions from z 2 -suppressed terms in the matrix element (i.e., twist-4 contributions). 

Contribution from the PW-term. The complete result for the Compton amplitude includes also the contribution 
from the PW term. Since this part of the matrix element of the contracted operator "decouples" and does not 
produce kinematical twist-3 terms, its contribution to the Compton amplitude should be transverse by itself. 



Indeed, a straightforward calculation, making use of the symmetry Eq.(4.7), gives 



1 D(gl) 

TU|pw-Term = 7 — r Vulu + %r v - 9^ (rq) + r^r v \ \ da- , (4-41) 

(rq) J_ x 1 - a 

which evidently satisfies 

<Z/j7)t;/|pW-Term = 0, TfjT^y | PW-Tcrm = 

0(t/q 2 ). (4.42) 

Hence this term can be regarded as a separate contribution. Alternatively, one may include the contribution from 
the PW-term into the basic SPD H(x, £) and all SPD's derived from it. Specifically, for £ > 0, the PW contribution 
to H(x,£) is D(x/C)0{\x\ < §; it contributes (£ - x)D'{x/£) 9{\x\ < £)/£ 2 [where D'(a) = (d/da)D(a)] to 
i?(x,£). Inserting these functions into Eqs. ( 4.29| ) and (1.52) one rederives Eq. ( }4.41 ). One can also observe that 



the PW term makes zero contribution to /(£), Eq. ( f4.35[) 
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D. Wandzura— Wilczek— type representation of twist— 3 SPD's 



We can express the results for the matrix elements and the Compton amplitude derived in Sections IV B and IV C 



in an alternative form, introducing new skewed distributions as Wandzura- Wilczek-type transformations of the 
original SPD's, H and A. This allows us to demonstrate that our approach is equivalent to the Wandzura- Wilczek 
approximation proposed in Ref. fl8[| . In this language, the singularity for transverse photon polarization described 
in Section IV C appears due to discontinuities of the new SPD's at x = £. 



Within our approach the step leading to a Wandzura- Wilczek-type formulation is to introduce in the parameter 
integrals over v the combination 



•0£ 



(4.43) 



as a new variable. With R(x, £) given by Eq.(4.30), we define a "transformed" function by way of 

J-i Jo 

R(v,t) 



dyR(y,£) dv5(yv + v£-x) 
x>0 /V M 



0(x < dy 

y-i J-i y-C 



where x ^ £. Note that the function 



R(x,0 



— [H(x,0-A(x,Q] 



(4.44) 



(4.45) 



is continuous at x = £. This important property comes about as follows: As can be seen from the graphical 
representation of the reduction formula, Fig. |l|, x = £ corresponds to the point where integration contour in 
Eq.(hl_4) passes through the upper corner of the "rhombus" representing the two-dimensional region of support 
of the double distribution, f(x,a). For a general double distribution, which not necessarily vanishes at a — ±\x\, 
the SPD's H and A individually will be continuous at x = £, while their first derivative in x will generally have a 
finite discontinuity. However, the difference, H — A, can be written as an integral of the type of Eq.( 4.14 ) but with 
a factor of 1 — a, which suppresses the contributions located in the corner of the integration region. As a result, 
not only the function H — A but also its first derivative in x are continuous at x = £. 

Matrix elements of vector-type string operators in terms of the transformed SPD R\y ■ In terms of the new 
function, Rw, we can express the matrix elements of the vector and axial vector operators, Eqs.( 4.1§| ) and (4.21), 
as 



(p-r/2\ [^(-z/2) 7 ^(z/2)] 



sym + total der. 



|p + r/2) 



dxe~ lx{pz) \ 2p, 



H(x,0-^A(x,0 



Act - p, 



(Az) 



(pz) 



[R w (x,0 - Rw(-x,£,)] 



(4.46) 



(p-r/2\ [^(~z/2) lal ^{z/2)} 



sym + total der. 



\p + r/2) 



^e aa0J -^rAf3p~t I dxe lx(pz) [R w (x, + Rw(~x, 0] ■ 



(4.47) 



Note that only the odd part of Rw(x,C) contributes to the vector operator, Eq.(4.46), and only the even part to 
the axial vector one, Eq.( 1.47 ). [We remind that the full vector matrix element includes also the contribution from 



the PW-term, Eq.( 4.22 ), which remains unchanged.] 

We can now establish the relation of our functions to the twist-3 SPD's introduced within the collinear expansion 
approach by Anikin et al. fl6|] . In Ref. fl6|| the matrix elements of the vector and axial vector operator were 
parametrized in the form (using our notation) 
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\ I ^e^(p-r/m-\n/2)\ ^ 
2 J 2lT { 7.75 



ip{Xn/2)\p + r/2) 



- e CTQ/ 3 7 n Q A/3P 7 H A {x,£). 



(4.48) 



where n is a light-like vector, n = 0, normalized as (pn) — 1, and (An) = 0. Setting in Eqs.(4T6|) and (4.47) 



z = An, so that the terms proportional (Az) drop out, and comparing with Eqs.( 4.48 ) we find that the Wandzura- 
Wilczek-type parts of these functions are given by 



H 3 (x,0 



w.w. 



(4.49) 
(4.50) 

(4.51) 

Compton amplitude in terms of the transformed SPD Rw ■ We can also express the Compton amplitude in terms 



= -j [Rw(x,0 - R w (-x,g)] 



of the new function, Eq.(4.44). This can be done either by computing the Compton amplitude as in Section IV C, 
using the parametrizations Eqs.( 4~46| ) and (4-47) for the matrix elements of the vector-type operators, or directly 
by rewriting the old result, Eq.( 4.29| ), in terms of the new SPD, Eq.( 4.44 ). We obtain 



T " V 2{pq) 



1 dx H ^ 
.! x-i + iO 



(9m + 3$J»m) a * 



■ X 



iO 



dx 

-1 



Rw(x,Q 



iO 



(4.52) 



To this one should add the contribution from the PW-term, Eq.(4.41), which remains unchanged. 

Singular integrals with the transformed SPD R\y- Let us now inspect the singular integrals in the new repre- 
sentation of the Compton amplitude, Eq.( 4.52j ). In particular, this will provide an alternative explanation of the 
divergence observed in Section IV C. In the second term, contributing to the amplitude for longitudinal photon 
polarization, the integral converges, since Rw{x,£) is continuous at x = — £. This follows immediately from the 
integral definition of Rw, Eq.(4.44), and the fact that the function R(x, £) in the integrand is bounded at x = — £. 



In the third term, however, contributing to the amplitude for transverse photon polarization, the integral with the 
singular function l/(x — £ + iO) would converge only if the function Ryy(x, £) were continuous at x = £. However, 
according to the definition, Eq.(4.44), Rw is given by the integral of the function R(y,£)/(y — £) from x to 1 if 
x > £, and from x to —1 if x < £. Hence there is no a priori reason why the limiting values of Rw{x,^) when 
approaching x = £ from below and from above should coincide. In fact, the difference of the two limiting values 
can be written as a principal value integral: 



flw(£ + o,0-iM£-o > 



1 . %Q 

dy r , 

-i y-€ 



(4.53) 



which is generally non-zero [in fact, this is just the real part of the integral Eq.( [t.35 )]. We thus conclude that in 
this term the x-integral is divergent. Thus, our conclusions regarding the existence of the singular integrals in the 
terms contributing for longitudinal and transverse polarization are the same as those made in Section IV C . This 
should be so; after all, the introduction of the Wandzura-Wilczek-type SPD, Rw, merely corresponds to a change 



of integration variables in the integrals in the second and third terms of Eq.( 4.29 ). Only the interpretation of the 
divergence is now different: It happens because of a discontinuity of the transformed function, Rw(x, £), at x = £. 
As we have seen, this discontinuity is built into the transformation at a very elementary level; it is not contingent 
on any dynamical properties of the underlying SPD's, H and A. This is in accordance with the fact that in the 



formulation of Section IV C the divergence in the u-integral is not influenced by the behavior of the SPD's. 
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V. CONCLUSIONS 



In this paper we have studied the DVCS amplitude at tree level, making use of the light-cone expansion in terms 
of QCD string operators in coordinate space. Let us summarize our main conclusions. 

First, we have found that transversality of the light-cone expansion can be maintained by including in addition 
to the twist-2 operators a minimal set of "kinematical" twist-3 operators, which appear as total derivatives of 
twist-2 operators. The "dynamical" twist-3 contribution from quark-gluon operators is not involved in restor- 
ing transversality of the twist-2 contribution; the contribution from quark-gluon operators should be separately 
transverse once also their total derivatives are included (this can be done by generalizing the iterative procedure 
described in Section [II B and Appendix [A]) . From the point of view of transversality alone it is thus a legitimate 
approximation to neglect this contribution.^ In this sense the original appeal of DVCS as a process describable 
in terms of a small set of twist-2 non-forward matrix elements, closely related to the parton distribution of DIS, is 
preserved. 

Second, we have found that retaining only the "kinematical" twist-3 contribution required by transversality 
leads to a finite result for the tensor amplitude corresponding to longitudinal polarization of the initial (virtual) 
photon. The tensor amplitude for transverse polarization contains a divergence already at tree level, whose origin 
are contributions with vanishing "longitudinal" component (in the parton model sense) of the momentum of the 
quark propagator. This divergence is specific to DVCS kinematics (technically speaking, they occur only if -q = £) 
and are absent for finite virtuality of the outgoing photon. However, we found that the divergent part of the 
tensor amplitude has zero projection on the polarization vector of the outgoing (real) photon. Thus, in both cases 
(longitudinally or transversely polarized initial photon) the physical hclicity amplitudes are finite. In this sense the 
approximation of keeping only the "kinematical" twist-3 contributions is consistent with the factorization of the 
DVCS amplitude. It may thus be used for estimating observables such as the DVCS cross section, the contribution 
from interference of DVCS with the Bethe-Heitler process, spin and charge asymmetries, etc. 

Our results can serve as a basis for computing kinematical power corrections to the DVCS amplitude, both 
of type t/q 2 and m 2 /q 2 , which are analogous to the well-known target mass corrections in DIS [^9 30 1. In the 
coordinate-space approach these corrections can be incorporated in a very elegant way, by including terms of order 
z 2 in the parametrizations of the matrix elements of the twist-2 operators. The latter are obtained by solving the 
harmonic condition on the matrix elements, Eq.( 3.15 ) [ fL3U3l| . 

In this paper we have computed the DVCS amplitude for the pion (more precisely, its isoscalar component), 
taking advantage of the fact that the pion has spin zero and definite C-parity. The generalization of the approach 
described in Section IV to the more realistic case of the nucleon is straightforward. The parametrization of the 
nucleon matrix elements of the twist-2 string operators in terms of double distributions and a PW-term has 
been discussed in Refs. In particular, this will allow to study the influence of the "kinematical" twist- 

3 contributions on the sum rules relating the nucleon SPD's in the forward limit to the quark orbital angular 
momentum in the nucleon Jl| . 
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sss'pjjjg situation with transversality in the light-cone expansion is very much reminiscent of the old problem of current 
conservation in the Bethe-Salpeter equation with reduced kernels, also known under the name of "exchange currents" fiiUiil . 
For a given truncation of the full Bethe-Salpeter kernel the "kinematical" modifications to the current operator, required 
to maintain transversality, can be computed from the truncated kernel in a well-defined way. In addition, however, there 
can appear "dynamical" exchange currents, which are separately transverse. See Ref. [|H) for a discussion of this problem. 
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of the "kinematical" twist-3 contributions is described. The divergences in the amplitude for transverse photon 
polarization were also noted recently by Kivel et al. [^0|. The generalization to the case of a nucleon target has 
been discussed in Ref. |5?j. 
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APPENDIX A: DECONSTRUCTING THE STRING OPERATORS 



In this appendix we demonstrate how the operator relations Eqs.(3.27) and ( 3.2S| ), expressing the vector string 
operator through its symmetric part and the total derivative of the string operator of opposite parity, can be solved 
recursively, giving rise to the closed-form expression stated in Eq.(3.43). As explained in Section IIIB, we shall 
not consider the quark-gluon operators appearing in the reduction process, since they are not needed to maintain 
transversality of the light-cone expansion. In principle, of course, generalizing the technique described below, the 
operator relations could be solved keeping also the quark-gluon operators. 

For brevity let us introduce the following notation for the basic (unsymmetrized) vector and axial vector string 
operators: 



Sa(z\X) 
S 5a {z\X) 



${X - z/2) 




ip(X + z/2). 



(Al) 



The R.H.S. of the relations Eqs.( 3.27 ) and ( 3.26 ) can be regarded as the action of a certain linear integral operator 
on the string operators: 



(CS) a (z\X) 



— 2 ^ a ctft-Y^cxRft 



dtt S~,(tz\X), 



(A2) 



where 



dX, 



(A3) 



With the help of this formal operator, Eqs.(3.27) and ( |3.28 ), with the quark-gluon operators neglected, can be 
written in compact form: 



(CS) a 
-(CS 5 ) C 



Q _ Q s .V m 

- s„ = S s J m . 



(A4) 
(A5) 



These equations decouple if one introduces instead of the vector and axial vector type string operators the left— 
and right-handed combinations (S a ± S , 5 CT )/2; the decoupling is, of course, a consequence of chiral invariance of 
QCD. The formal solution of the system Eqs.(A4) and JAE| ) is 

1 



S§o- — 



i-c 2 
1 

T^c 2 



gsym 

rtsym 
Ok 



1-C 2 
C 

1-C 2 



gsyn 



(A6) 
(A7) 



These relations allow to express the full string operators in terms of totally symmetrized (i.e., twist-2) operators. 
The inverse of the linear operator 1 — C 2 here is understood in the sense of a series expansion: 



1 



1-C 2 

c 

1-C 2 



= l + £ 2 - 
= C + C 3 



£ 4 + 
-C 5 + 



(A8) 
(A9) 



By iterating the operator C one can easily convince oneself that the even and odd powers of C observe the following 
patterns (n > 1): 



(C 2n S) c 
(£ 2n - 1 5) c 



(-1) 



n-l 



2(2n- 1)! 
2(2n-2)l' 



[R„z 1 - (Rz)g ai \ 

xfj-fZaKp 



(Rz) 



2 

2n-2 



dv vv 2 "- 1 



SJvz), 



dv 



Sj(vz). 



(AlO) 
(All) 
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These expressions apply up to terms involving operators of the form 



(A12) 



whose matrix elements in DVCS kinematics are proportional to t or m 2 . More rigorously, Eqs.( A10| ) and (All) 
can be proven be induction. Obviously, Eq.( All ) with n = 1 describes the first power of the operator, C, Eq.([A2|). 
Acting on Eq.(All) (for general n) with C we have to remember that Eq.(A2) implies that the distance, z, is to be 
rescaled by the parameter, t, everywhere, not only in the arguments of the fields, so that we pick up a total factor 
oft 2 "" 1 . We thus obtain: 



[C{C**- x S)] a = (C 2n S) a = 



4(2n-2)! 



(Rz) 



2n-2 



dtt 



2n 



RpRp 



duuu S^/(tuz). 



Using the well-known expression for the contraction of two e-tensors we find 

^crap-y^^a' P'-i' z a z a' RpRfi' = [R<j z ~f' ~ (R z )9cr~f'} terms OC R 2 , R^i . 



(A13) 
(A14) 



The parameter integral in Eq.(Al3) can be simplified by introducing a new variable, v = tu: 



dtt 2 



dt 



duu(\-u) 2n - 2 Sy(tuz) 



dvv(t-v) 2n - 2 S Y (vz) 



dv / dtv(t-v) 2n - 2 Sy(vz) 



1 



2n- 1 



dvvil-v) 2 "- 1 S Y (vz). 



(A15) 



Combining these statements we see that Eq.( A13| ) is indeed equal to the proposed expression, Eq.( A10| ). In an 
analogous way one shows that acting on Eq. ( A1C ) with C one reproduces Eq. ( All ) with n replaced by n + 1 . 

From Eqs.( A10) and (All) we see that the series of even and odd powers of the operator, C, Eqs.(|A8|) and (|A^), 
have the form of the power series of the sine and cosine functions, respectively: 

2n-l 



(-1)' 



^ (2,/ I )' 

n— 1 v ' 

g (_l)«-i 



(Rz) 



(Rz) 



.(Rz) 



2n~2 



— COS 



2 

(Rz) 



^ (2n-2)! 

n— 1 v ' 

Substituting this in Eqs.(|A6|), (^) we finally obtain [cf. Eq.( |3.43[ )]: 

1 f 1 

S a (z) = S s J m {z) + -[RaZ-y - {Rz)g a7 } / dvv sin 

2 Jo 

(Rz) 



(A16) 
(A17) 



(Rz) 



s s y m (vz) 



2 £(7a{3yZ a Rfi 



dv v cos 



Sl* m (vz) 



(A18) 



and an analogous result with S and S$ interchanged. Note that this expression applies up to operators of the form 
Eq.(A12), which give power-suppressed contributions to the DVCS amplitude. 

We remind again that the full expression for the string operator contains also quark-gluon string operators, 
which upon symmetrization produce operators of twist > 3. We have not kept track of them, since they turn out 
not be essential for restoring transversality of the twist-2 contribution. In principle they could be included in the 
above iterative procedure. 
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